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Abstract 


Current  research  in  partial  differential  equatioris  is  extensive, 
varied  and  deep.  A  single  lecture,  if  it  is  not  to  be  a  mere 
catalogue,  can  present  only  a  partial  list  of  recent  achievements,  some 
comments  on  the  modern  style,  i.e.  the  kinds  of  problems  chosen  and 
methods  used  for  solution,  and  cautious  speculations  on  future  trends. 
The  choice  of  examples  is  of  course  shaped  by  the  personal  taste  of  the 
speaker  and  limited  by  his  expertise. 

The  first  part  of  this  lecture  is  such  an  overview;  it  is  followed 
by  a  more  detailed  discussion  of  two  topics  with  which  the  speaker  has 
some  familiarity,  one  concerning  a  linear,  the  other  a  nonlinear 
problem  in  partial  differential  equations. 


XV 


la)  Linear  Problems.  In  the  last  few  years  a  number  of  problems 
concerning  linear  partial  differential  operators  on  manifolds  with 
boundaries  have  been  solved  or  are  nearing  solution.  Thanks  to  the 
researches  of  Melrose  [29],  Taylor  [36],  Ivrii  and  others  we  understand 
well  the  propagation  of  signals  along  reflected  glancing  and  gliding 
rays,  the  clue  to  many  problems  in  diffraction  and  scattering. 
Microlocal  analysis,  the  modern  version  of  wave-ray  duality,  has 
provided  the  tools:  pseudo-differential  operators,  Fourier  integral 
operators,  Harailtonian  flows  and  Lagrange  manifolds.  In  his  recent 
work  Charles  Fefferraan  [14]  makes  use  of  a  sophisticated  version  of  the 
uncertainty  principle.  Another  versatile  modern  technique  is  the  use 
of  trace  formulas  to  link  spectral  and  geometric  information. 

For  a  thorough  documentation  of  the  sucesses  of  the  modern  theory 
of  linear  partial  differential  equations  we  have  to  await  the 
publication  of  Iformander's  3-volurae  treatise  [18],  but  it  is  clear  that 
the  successes  have  been  so  sweeping  that  they  have  radically  altered 
the  course  of  research  in  this  field.  I  believe  that  in  the  future  we 
shall  see  more  applications  of  the  methods  and  results  of  the  theory  of 
linear  partial  differential  equations  to  other  fields  of  mathematics; 
examples  from  the  past  are  of  use  of  PDE  methods  in  several  complex 
variables  and  quasiconf orraal  mappings.  We  are  likely  to  see  more 
special  questions  raised,  from  sources  inside  and  outside  mathematics, 
and  more  detailed  answers  given;  mere  preoccupation  with  existence  and 
uniqueness  questions  is  likely  to  diminish. 

Section   2   contains   a   brief   description  of  wave  propagation  on 

complete  manifolds  of  constant  negative  curvature.   One   of   the   tools 

,  used   Is   the   Radon   transform,   less   popular  than  its  more  glamorous 


sister,  the  Fourier  transform,  but  more  appropriate  in  some  situations. 
See  also  [17], 

b)  Nonlinear  probleris. 

The  strides  that  have  been  made  recently  in  the  theory  of 
nonlinear  PDE-s  are  as  great  as  in  the  linear  theory.  Unlike  the 
linear  case,  no  wholesale  liquidation  of  broad  classes  of  problems  has 
taken  place;  rather  it  is  steady  progress  on  old  fronts  and  on  some  new 
ones,  the  complete  solution  of  some  special  problems,  and  the  discovery 
of  some  brand  new  phenomena.  The  old  tools  —  variational  methods, 
fixed  point  theorems,  degree  of  mapping  and  other  topological  methods 
—  have  been  auguraented  by  some  new  ones.  Preeminent  among  these, 
indi^spensable  for  discovering  new  phenomena,  is  numerical 
experimentation;  but  it  is  likely  that  in  the  future  numerical 
calculations  will  be  part  of  proofs. 

We  shall  discuss,  very  briefly,  three  topics: 

i)  Viscous,  incompressible  flows. 

ii)  Hyperbolic  systems  of  conservation  laws  and  shock  waves. 

iii)  Completely  integrable  systems, 

i)  Viscous ,  incompressible  flows 

In  spite  of  a  claim  by  Kaniel,  [20],  laid  to  rest  by  D.  Michelson, 
the  existence  for  all  time  of  strong  solutions  of  the  Navier-Stokes 
equation,  and  the  uniqueness  of  weak  solutions,  in  three  dimensional 
space  are  very  much  open  questions.  Ue  have  learned  more  about  the 
singularities  of  weak  solutions,  in  particular  about  the  Hausdorff 
dimension  of  the  singular  set.  Already  Leray  has  shovm  that  every 
solution   Is   continuous   If   we   eliminate  a  closed  set  of  t  with  zero 


Rausdorff  measure  of  dimension  1/2.  B.  Mandelbrot  has  raised  the 
question  of  what  the  Hausdorff  dirneasion  of  the  possible  singularities 
of  weak  solutions  is  in  space  and  titne.  The  first  results  on  this 
important  question  were  obtained  by  V.  Scheffer  [33];  the  latest  word 
is  the  following  theorem  of  Caffarellis  Kohn  and  Nirenberg  [5]: 

The  one  dimensional  Hausdorff  measure  of  the  set  of   singularities 
of  a  suitable  weak  solution  in  x,t-space  is  zero. 

Turbulence,  surely  one  of  the  outstanding  problems  of  mathematics, 
can  be  described  by  the  long-time  behavior  of  typical  solutions  of  the 
NS  equations.  When  viscosity  is  large  compared  to  the  force  driving 
the  flow  there  exists  exactly  one  stationary  flow  to  which  all  flows 
tend.  As  the  force  is  increased,  this  stationary  flow  becomes 
unstable,  i.e.  any  slight  perturbation  drives  it  to  another,  stable, 
stationary  flow.  When  the  force  is  increased  still  further,  this  too 
becomes  unstable  and  the  flow  tends  to  yet  another  stationary  flow  or 
possibly  to  a  periodic  flow.  As  the  force  is  increased  further  the 
flow  becomes  more  and  more  chaotic.  This  chaotic  flow  is  concentrated 
around  a  so-called  attractor  set,  i.e.  a  set  consisting  of  points  of 
accumulation  of  a  single  flow  driven  by  a  force  that  is  independent  of 
t.  Such  sets  are  invariant  under  the  Navier-Stokes  flow;  concerning 
these  Foias  and  Temam  have  proved  the  following,  see  [16]: 

A  bounded  set  that  is  invariant   under   the   strong  Navier-Stokes 
flow  in  a  bounded  domain  has  finite  Hausdorff  dimension. 

The   dimension   of   such   sets  may  go  to  infinity  as  the  viscosity 
tends  to  zero. 

Further  results   along   these   lines   have   been   obtained   by   P. 
Constantin  and  C.  Foias. 


The  simplest  testing  ground  for  ideas  of  instability  and 
turbulence  of  viscous  fluids  are  the  Couette-Taylor  flows,  i.e.  flows 
between  two  concentric  cylinders,  the  inner  one  rotating  with  some 
angular  velocity  oj .  If  the  cylinders  have  infinite  length,  then  there 
is  a  stationary  flow  that  is  independent  of  the  angle  0  and  distance  z 
along  the  axis  of  the  cylinders.  For  id  low  enough  this  flow  is  stable; 
as  CO  increases,  this  flow  becomes  unstable,  yielding  stability  to 
another,  z-dependent,  flow  consisting  of  a  stack  of  Taylor  vortices, 
named  after  their  discover.  As  lo  increases  further  this  flow,  too, 
becomes  unstable  and  gives  way  to  a  0-dependent,  periodic  flow, 
further  increase  in  w  leads  to  more  and  more  complicated  flows. 

There  is  a  wealth  of  experiiiental  studies  of  Couette-Taylor  flows, 
revealing  a  bewildering  variety  of  steady  and  unsteady  flows,  see  e.g. 
Benjamin,  [2].  The  understanding  of  these  (which  must  also  take  into 
account  the  finiteness  of  the  cyclinders)  is  a  profound  challenge  to 
theoretical  and  computational  fluid  dynamicists. 

Flows  without  any  driving  force  to  maintain  them  decay  because  the 
viscous  forces  dissipate  energy.  Recently  Foias  and  Saut  [15]  have 
shown  that  the  rate  of  decay  is  exponential,  the  same  as  for  the 
corresponding  linearized  Stokes  flow;  they  have  further  shown  that  the 
Stokes  and  the  Navier-Stokes  flows  are  linked  by  a  wave  operator. 

When  viscosity  is  zero,  as  in  the  Euler  equation,  no  imposed  force 
is  needed  to  maintain  the  flow.  Existence  and  uniqueness  is  known  in  2 
dimensions  but  is  doubtful  when  n  =  3.  Extensive  calculations  by  S. 
Orszag  and  his  colloborators  on  the  Taylor-Green  vortex  problem,  [30] 
reveal  a  bewilderingly  complicated  flow;  as  time  goes  on,  smaller  and 
smaller  scale  features  appear  until  the  numerical  method  —  a   spectral 


method  keeping  track  of  more  than  100  million  Fourier  coefficients  — 
is  unable  to  resolve  thera.  Another  set  of  calculations  by  Orszag 
employs  the  Taylor  series  in  time,  up  to  order  88,  summed  in  a  cunning 
fashion  to  elude  singularities  in  the  complex  t-plane;  the  features 
revealed  in  the  two  calculations  are  similar. 

Another  set  of  calculations,  not  nearly  so  machine-intensive  as 
Orszag's  has  been  performed  by  Chorin,  see  e.g.  [7].  He  considers  an 
initial  value  problem,  periodic  in  space,  where  vorticity  is  initially 
confined  to  a  narrow,  slightly  crooked  tube.  The  basic  variable  is 
vorticity,  and  the  calculation  takes  into  account  that  the  vorticity  is 
confined  to  the  tube,  which  stretches  and  twists  with  the  flow.  Using 
a  number  of  bold  simplifications  the  calculation  is  carried  out  long 
enough  to  indicate  that  after  a  finite  time  the  vortex  tube  will  be 
stretched  so  thin  that  its  Hausdorff  dimer^sion  becomes  "2.5,  a 
prediction  of  Mandelbrot's,  [28].  Another  calculation  by  Chorin, 
employing  a  rescaling  reminiscent  of  the  renormalization  group  of 
physicists,  leads  to  the  same  conclusion. 

ii)  Hyperbolic  systems  of  conservation  laws  and  shock  waves. 

A  system  of  n  conservation  laws  in  one  space  variable  x  and  in  t 
is  of  the  form 

"t  +  f(^)x  =  0. 

u(x,t)  in  E^;  the  system  is  strictly  hyperbolic  if  the  matrix  Vf(u)  has 
real  eigenvalues  for  every  u  in  E'^. 

The  basic  problem  is  the  initial  value  problem:  given  u(x,0)  = 
'Jq(v)j   show  the  existence  of  a  solution  u(x,t)  for  all  t,  in  the  class 


of  discontinuous  functions,  satisfying  the  conservation  law  in  the 
sense  of  distributions,  and  an  entropy  condition,  of  the  form 

s^  +  gj^  <  0, 

where  s  =  s(u)  is  an  entropy,  g  =  g(u)  entropy  flux.   The  two  satisfy 

VsVf  =  Vg  , 

and  s  is  required  to  be  a  convex  function  of  u. 

Numerical  evidence  indicates  strongly  that  various  difference 
schemes  for  solving  conservation  laws  converge;  yet  until  recently  no 
proof  had  been  given  for  systems  .with  more  than  one  state  variable  u. 
Similarly,  physics  strongly  suggests  but  mathematics  had  been  unable  to 
prove  that  if  u  =  u  (x,t)  solves  the  viscous  equation 

u^.  +  f(u)^=  e  D(u)xx  .  e  >  0. 

u(x,0)  =  Uq(x), 

D  an  appropriate  n  x  n  viscosity  matrix,  then  as  e  tends  to  0,  u^(x,t) 
converges  to  a  solution  of  the  system  of  conservation  laws  that 
satisfies  an  entropy  condition. 

This  year  Ron  Di  Perna  [11],  succeeded  in  proving  such  convergence 
theorems  for  the  equation  governing  the  flow  of  an  isentropic  gas 
satisfying  a  polytropic  equation  of  state,  with  the  artificial 
viscosity  D  H  I.  Among  the  many  ingredients  are  two   beautiful   general 


ideas   of   Tartar  and  Murat ,  [35].   One  Is  a  characterization  of  strong 
convergence  In  terms  of  weak  convergence: 

Suppose  u.(v)  is  a  sequence  of  mappings  from  t  to  ^^',  uniformly 
bounded  in  L" ;  an  then  there  is  a  subsequence  such  that  for  every 
continuous  function  g  in  1™  tlie  weak  limits 


w-llm  g(u.(y)) 


exists.      These    lira.lts    can   be    represented   as 


/g(u)    dVy(u)    , 


where  v  is  a  probability  measure  in  S.™  parametrized  by  y  in  E.™.  The 
subsequence  u.  converges  strongly  if  the  measures  v  have,  for  each  y, 
a  single  point  as  support. 

The  second  idea  is  compensated  compactness:  Let  v.  and  w-  be  two 
sequences  of  mappings  from  V^  to  iR  ,  if  both  converge  weakly  in  the  L2 
sense  to  v  and  w  respectively,  and  if  div  v,  and  curl  w-i  lie  in  compact 
sets  in  the  ^l~]_^j^  topology,  then 


lim  J  v.  •  w-j  dy  =  /  v  •  w  dy. 


Tartar  himself  has  used  these  ideas  to  give  new  proofs  of  known 
convergence  tlieorems  for  scalar  conservation  laws;  Di  Perna  has  shown 
how  to  use  them  for  systems  with  two  variables. 

Very  little  is  known  about  existence  of  discontinuous  solutions  in 
more  than  one  space  variable;  even  short  time  results  are  of  recent 
origin,    see    A.   Majda's   talk   at   this   Congctjss.    Yet   numerical 


calculations,  done  with  care  and  ingenuity,  see  e.g.  Colella  and 
Woodward,  [9],  converge  and  give  solutions  consistent  with 
experieraents. 

In  one  space  dimension  we  know,  at  least  for  simple  systems,  that 
shock  formation  and  interaction  severely  limit  the  amount  of 
information  that  a  one  dimensional  flow  field  can  contain.  Something 
similar  mist  be  true  in  higher  dimension,  but  the  mechanism  causing  it 
is  not  understood. 

In  his  report  to  this  Congress,  S.  Klainerman  will  describe  some 
recent  results  on  long  term  existence  of  regular  solutions  of  the 
Initial  value  problem  for  non-linear  hyperbolic  equations  in  several 
space  variables. 

The  question  of  uniqueness,  subject  to  an  entropy  condition,  is 
not  satisfactorily  settled  even  in  one-space  dimension,  not  even  for 
the  equation  of  compressible  flow,  in  spite  of  the  important  pioneering 
work  of  Oleinik,  and  more  recent  work  of  Di  Perna. 

There  are  intriguing  open  problems  concerning  ,  transonic  flows 
around  given  contours,  with  volocity  prescribed  at  infinity.  An 
ingenious  method  of  Garabedian  yields  large  families  of 
aero-dynamically  interesting  smooth  flows,  but  a  basic  theorem  of 
Morowetz  shows  that  for  a  generic  contour  no  shockless  flow  exists. 
The  basic  problem  is  to  prove  the  existence  of  a  flow,  with  shocks,  and 
to  prove  its  uniqueness.  Recent  numerical  calculations  of  A.  Jameson 
indicate  that  in  the  potential  flow  of  approximation  there  may  be  many 
solutions. 

iii)  Completely  integrable  systems 

This  chapter  In  mathematics,  barely  15   years   old,   continues   to 


fascinate  analysts  and  algebraists,  as  well  as  physicists.  The  effort 
has  been  truly  international,  and  has  paid  off  in  the  discovery  of  new 
completely  integrable  systems,  many  of  physical  interst,  some 
containing  two  space  variables.  The  algebraic  classification  of  these 
systems  has  progressed,  and  new  connections  with  other  branches  of 
mathematics  and  physics  have  been  found,  such  as  the  t -function  of 
Sato,  Mwa  and  Jimbo,  [32],  Three  books  have  appeared  recently  on 
solitons  and  scattering  theory  [1],  [6]  and  [40],  and  the  work  of  Beals 
and  Coifman,  [3].  The  speaker  will  restrict  his  report  to  a  few 
scattered  remarks  on  the  analytic  side  of  the  matter. 

a)  Solutions  of  completely  integrable  partial  differential 
equations  lie  on  infinite  dimensional  tori.  Numerical  experiments  with 
such  equations,  see  e.g.  [21],  furnish  numerical  approximations  that 
appear  to  lie  on  tori,  necessarily  finite  dimensional.  THis  indicates 
that  some  infinite  dimensional  analogue  of  the  KAJi  theory  might  be 
true;  no  such  result  is  known. 

b)  The  sine-Gordon  equation 


"tt  -  u^jj  +  sin  u  -  0    ^  U 


has     explicit  solutions  that  decay  exponentially  in  x  as  ]x|  -»■  =>  and 
(^r',;  periodic  in  time.   If  the  function  sin  u  is  replaced  by  g(u),   Coron 

[10]   has   shown  that  no  such  solution  can  exist  when  the  time  period  T 

211 
is  <   ; .   •  Coron  and  Brezis  conjecture   that   there   are   no   periodic 
g  CO) 

solutions  of  any  period,  except  for  very  special  functions  g. 

c)   The  explicit  solution  of  the  initial  value  problem  for  the  KdV 
equation 


"t  -Su^x  +  ^'^   "xxx  =  0 


u(x,0)  =  u(x) 


in  terms  of  the  scattering  transform  makes  it  possible  to  determine 
explicitely  the  limit  of  the  solution  u(x,t,E)  as  e  >  0.  This  rather 
interesting  limit  is  described  in  Section  3. 

These  altogether  too  brief  remarks  on  nonlinear  PDE  were  confined 
mostly  to  problems  arising  in  mathematical  physics;  it  is  the  richest 
source  of  such  material,  but  not  the  only  one:  geometry  is  another,  see 
S.T.  Yau's  report  to  the  Helsinki  Congress  [38].  The  speaker  has 
neither  the  knowledge  nor  the  time  to  report  on  progress  in  this  very 
active  area  in  the  last  five  years,  but  he  cannot  resist  mentioning  the 
very  recent  demonstration  by  Brezis  and  Coron  of  the  existence  of  two 
surfaces  of  constant  mean  curvature  spanning  a  prescribed  plane  curve, 
not  too  large;  the  proof  is  a  marvel  of  subtlety. 

There  hasn't  even  been  enough  time  to  mention  all  the  subjects  in 
mathematical  physics  that  have  been  traditionally,  but  especially  in 
the  recent  past,  rich  sources  of  problems  in  nonlinear  PDE:  elasticity 
theory,  electromagnetic  theory  and,  more  recently, 
magnetohydrodynamics.  Two  topics  which  need  more  help  from 
mathematicians  than  they  are  getting  now  are  multiphase  flow  and 
combustion  see  e.g.  [8]  and  [27].  In  both  it  is  of  great  importance 
to  understand  the  nature  of  turbulent  regimes;  but  in  multiphase  flow, 
as  in  aero  and  hydrodynamics,  turbulence  is  detrimental;  in  combustion 
it  is  beneficial  since  It  promotes  the  mixing  of  fuel  and  oxydi^.er.  On 
the  other  hand  Shockwaves  are  detrimental  for  combustion,  since  they 
produce  entropy  which  decreases  the  efficiency,  of  conversion  of  heat 
into  mechanical  energy. 

10 


2.  The  Laplacc-Beltrami  operator  on  complete  RieTnanulan  manifolds   with 
constant  negative  sectional  curvature. 

In  a  series  of  papers,  [23]  -  [26],  R.  S.  Phillips  and  the  speaker 
have  analysed,  fairly  completely,  the  spectral  properties  of  the 
Laplace-Be Itranii  operator  on  manifolds  F  as  above  in  the  case  when  F 
has  infinite  volume  and  is  geometrically  finite.  This  extends  to  all 
dimensions  the  previous  work  of  Patterson,  [31],  in  the  case  n  =  2,  and 
allows  F  to  have  cusps  of  all  kinds. 

The  universal  cover  of  F  is  hyperbolic  space  B  ;  F  itself  can  be 
identified  with  the  quotient  ^^/T  ,  T  a  discrete  subgoup  of  isoraetries 
of  E^,  More  concretely  F  can  be  identified  with  a  fundamental 
polyhedron  F  of  E^  mod  T,  Conversely,  if  F  is  a  discrete  subgroup  of 
the  group  of  all  isoraetries,  then  l^/F  =  F  is  a  complete  Riemannian 
manifold  with  constant  negative  sectorial  curvature;  if  T  contains 
elliptic  elements,  F  has  harmless  singularities  along  submaaif olds. 

We  shall  use  the  Poincare  model  for  I^,  i.e.  the  upper  half  space 
(x,y),  X  in  K"^"^,  y  >  0,  equipped  with  the  metric 


ds2  =^]^L±^ 

y2 


The  set  of  points  at  infinity,  (g  ,0)    » ,  is  denoted  by  B. 

The  L^   norm  and  Dirichlet  integral,  invariant  under  isometries  of 
Sy,  will  be  denoted  by  H(u)  and  D(u): 


H(u)  =/  !u2|  1^,    D(u)  =/(|uj2  4-  lu  |2)l!i^ 
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The    invariant  Laplace-Be  It raral   operator  Lg   is   defined   in   terras    of    these 
quadratic   forms: 

H(u,Lqv)    =  -D(u,v) 
for  all  C^   functions   u   and   v, 


Lq  =  y2    (A49y2)    -    (n-2)y  3y    , 


where  A  =  \  dx  is  the  Euclidean  Laplace  operator.  Using  the 
Friedrichs  extension,  Lg  can  be  made  into  a  nonpositive  self  adjoint 
operator  with  respect  to  H. 

Similarly  we  denote  the  L'^  and  Dirichlet  integrals  over  F  by  Hp(u) 
and  Dp(u);  here  u  is  any  Cg  autonorphic  function  with  respect  to  a 
given  discrete  subgroup  T  of  isometries,  and  F  is  a  fundamental 
polyhedron  for  F.  In  what  follows  we  assume  that  F  has  a  finite  number 
of  sides,  i.e.   that  T    is  geometrically  finite. 

Discrete   subgroups   can  be  classified  by  the  geometric  properties 
of  their  fundamental  polyhedra  F  into  the  following  classes: 
i)  F  is  compact 
ii)  F  is  noncompact  but  has  finite  volume: 


V(F)  =/p£^< 


iii)    F    has    infinite   volume. 
The  spectral    properties    of      Lg     are      sharply      different      in      these 
cases. 
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In  case  i)  it  follows  by  standard  elliptic  theory  that  the 
spectrum  of  Lq  ±s  standard  discrete,  i.e.  pure  point  spectrum 
accumulating  only  at  -« .  The  present  work  is  concerned  nainly  with 
case  iii) ;  our  results  are: 

a)Lg  has  absolutely  continuous  spectrum  of  infinite  multiplicity 
in  (-<^,  -(Illi)2. 


b)   Lg   has   at   most   a   finite   number  of  point  eigenvalues,  all 
located  in  the  interval  (-("~  )^,0). 


c)  Lg  has  no  singular  spectrum;  even  the  point  spectrum  may  be 
empty.  However,  Beardon  and  Sullivan  [34]  have  shown  that  if  F 
contains  a  cusp  of  highest  rank,  then  there  is  at  least  one  point 
eigenvalue.   We  have  a  new  proof  of  this  result.  Theorem  6.4  in  [25]. 

Jorgensen, [1 9 ]  has  constructed  interesting  examples  of  groups  of 
isometries  of  B^  whose  fundamental  polyhedron  has  infinitely  many 
sides.  For  these,  Epstein,  [12]  has  shown  that  Lg  has  infinite 
dimensional  spectrum  in  (-1,0). 

Case  ii)  is  a  curious  mixture  of  i)  and  iii):  Lg  has  absolutely 
continuous  spectrum  in  (-«> ,-(  )  ) ,  but  only  of  finite  multiplicity, 
which  is  equal  to  the  number  of  cusps  .  There  is  no  singular  spectrum 
but  there  is  a  point  spectrum,  accumulating  only  at  ^°.  In  many 
special  cases  it  is  known  that  this  point  spectrum  is  ample.  In  the 
general  case,  Selberg  has  established  a  relation  between  the  density  of 
the  point  spectrum  and  the  winding  number  of  the  determj.riant  of  the 
scattering  matrix;  see  also  page  205-216  of  [23].   To  give  an  absolute 
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estimate  of  the  number  of  point  eigenvalues  remains  a  challenging  open 
problem. 

We  return  now  to  case  iii).  Earlier  s'tudies  of  the  continuous 
spectrum  of  Lq  proceeded  by  constructing  explicitly  a  spectral 
representation  of  Lg;  the  generalized  eigenf unctions  of  Lg  entering 
this  spectral  representation  are  Elsenstein  series,  constructed  by 
analytic  continuation. 

Our  approach  is  entirely  different;  it  is  applicable  to 
representing  operators  whose  continous  spectrum  has  uniform 
multiplicity  on  the  whole  line.  Let  A  be  an  anti-self  adjoint  operator 
whose  spectrum  is  of  uniform  multiplicity  on  the  whole  imaginary  axis. 
Then  the  spectral  representation  for  A  can  be  thought  of  as 
representing  the  underlying  Kilbert  space  H  by  L~(t,N),  N  some 
auxiliary  Hilbert  space  whose  dimension  equals  the  multiplicity  of  the 
spectrum  of  A; 

Each  f  in  H  is  represented  by  an  h^  function  K(a),  o  in  1,  the 
values  of  K  lying  in  N.  Since  A  is  antiself ad joint,  Af  is  represented 
by  iaK(a). 

Denote  U(t)  the  unitary  group  whose  generator  is  A:  U(t)=exp  tA. 
Then 

U(t)f  ^>  e^o^KCa)  . 

The  Fourier  transform  of  this  representation  with  respect  to  a 
gives  another  representation  of  H  by  L'^(1.,N),  where  each  f  in  H  is 
represented  by 

f  -►->  k(s)  ,  k  =  K  . 
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Then 

A  f  -H-^  -A  k(s) 
as 

and 

U(t)  f  ■<-+  k(s-t)  ,  U(t)  =  exp  tA. 

This  is  called  translation  representation. 

Of  course,  conversely,  the  Fourier  transform  of  a  translation 
representation  is  a  spectral  representation. 

We  show  now  how  to  construct  a  translation  representation  for  the 
unitary  group  associated  with  the  non-Euclidean  wave  equation 

"tt  -  ^"  =  °  ' 
where 

Note  that  if  a),  b) ,  c)  hold  then,  apart  of  the  finite  point  spectrum, 
L  has  continuous  spectrum  of  uniform  multiplicity  on  (-o°,0). 

The  group  associated  with  the  wave  equation  consists  of  the 
operators  mapping  initial  data  into  data  at  time  t: 


U(t)  :  {u(0),  Uj.(0)}  >  {u(t),  ut(t)} 


The  generator  of  U  is 


^0  L) 
A{u,Uj.}  =  {u^.u^t-J  =  Wf^  "^  =  tu,ut.}5i  ol 
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0  TO 

Note  that  A-  =  (  ^^^  ^^^^  ^  having  continuous  spectrum  of  uniform 
multiplicity  on  1_  is  consistent  with  A  having  continuous  spectrum  of 
uniform  nultiplicity  on  i  5.. 

Most  properties   of   the  non-Euclidean  wave  equation  follow  from 

standard  hyperbolic  theory: 

i)  The  initial  value  problem  is  properly  posed. 

ii)  Signals  propagate  with  speed  <  1 

iii)  If  the  initial  data  are  automorphic,  so  is  the   solution  for 
all  t. 

iv)  Energy  is  conserved,  where 


E  =  H(u^)  -  H(u,L  u)  =  H(uj.)  +  D(u)  -  (^)^  H(^)  • 


Finally  we  have  the  special  property 

v)  For  n  i^vj  J  the  Huygens  property  holds,  i.e.  signals  propagate 
with  speed  =  1 . 

It  is  not  hard  to  show  by  integration  by  parts  that  for  Cq°°  data 
in  ffl  ,  the  energy  E  is  positive.  We  denote  by  H  the  completion  in  the 
energy  norm  of  Cq"  data.  It  follows  from  conservation  of  energy  that 
U(t)  is  unitary  for  the  energy  norm. 

We  define  the  Radon  transform  of  a  function  u  in  H^  by 

a  -     i  u  dS. 

5(s,B) 

Here  5  (s  ,B  )   is   the   horosphere   centered   at  the  point  B  at  " ,  whose 
distance  from  the  origin  is  s.  It  is  well  known  that 
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Lu=e  Illis  a2  e^s   a   . 


Now   take   the  Radon   transform  of   the   wave   equation: 


0  =  a^t  =  La  -  a^t  -  e—r^  dl  -r^  a 


n-l_      .    n-1 
s 


Introduce 


n-1 
Q  =  -y*   V    ; 

then 

°   =  ^tt   -  ^ss   =   Ot  +  9s)    (^t   -  ^s). 

from  which   it    follows   that   v^   -  Vg ■ is   a   function  of   s   -     t.      We      define 
now 

R+{u,Ut.}    =  P(vt   -  Vg)    =  P(e^   a^   -  Sge^Q), 

P  an  appropriately  chosen  operator  in  s  that  commutes  with  translation. 

R+  is  a  translation  representation  of  H  for  V(t),  i.e. 

i)  R+  {u(t),  Ut.(t)}  (s)  =  R+  {u(0),  Uj-(O)}  (s-t) 

ii)  /(R+{u,Ut.})2  ds  dm(6)  =  E(u) 

iii)  R^  maps  H  onto  L2  (E  •  B). 

Of  course  i)  follows  from  the  way  R^  was  constructed;  for  ii)  and 
iii),  see  Theorem  3.4  [24]  when  n  =  3  and  P  =  c9  . 

We  turn  now  to  the  automorphic  case.  Here  energy,  defined  as 
before  by 
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Ep(u)  =  Hp(ut)  +  Dp(u)  -  (I^)\(u)  ,. 


Is  not  necessarily  positive.  It  was  shown  in  Section  4  of  [25]  that 
one  can  add  a  quadratic  function  K(u)  to  the  energy  so  that  G(u)  = 
Ep(u)  +  K(u)  is  positive,  and  that  K  is  compact  with  respect  to  G.  It 
follows  from  this  that  if  Ep  is  negative  on  a  subspace,  that  subspace 
is  finite  dimensional.  Since  for  u.  =  0  the  energy  is  -(u,L  u  )p,  it 
follows  that  the  positive  spectrum  of  L  consists  of  a  finite  number  of 
eigenvalues.  It  can  be  further  shown,  using  the  fact  that  F  contains 
full  neighborhoods  of  points  at  »,  that  L  has  no  negative  eigenvalues, 
see  Theorem  4^8  of  [25];  this  is  a  non-Euclidean  version  of  a  classical 
result  of  Rellich  and  Vekua. 

It  follows  from  the  form  Ep  =  Hp(uj^)  -  Hp(u,L  u)  that  if  u  is 
orthogonal  to  all  the  eigenf unctions  of  L,  then  Ep  >  0.  If  u(0)  and 
u  (0)  are  both  orthogonal  to  all  the  elgenfunctlon  it  follows  that  so 
is  u(t),  for  any  t.  We  denote  this  space  of  initial  data  by  H  . 
Clearly  H  is  invariant  under  the  solution  operators  Up(t)  for  the 
automorphic  solutions  of  the  wave  equation. 

We   define  now  a  translation  representation  Rjj.  of  H   for  the  group 

Up(t);  it  has  M  +  1  components,  M  being  the  number  of  cusps  of   maximal 

p 
rank.    The   zeroth   component  of  R^  is  R.,.,  defined  as  before;  each  of 

the  rem,ainlng  components  are  associated  with  cusps  of  maximal   rank   as 

follows: 

Map   the   cusp  to  the  point  " ,  so  that  it  has  the  form  E^  x  (a,"), 

Eoo  a   fundamental   polyliedron   In  Euclidean   space   of   the   parabolic 

subgroup   keeping  ■^   fixed.   Since  the  cusp  Is  of  ir.iiximal  rank,  F^  has 

finite  volume.   Denote  by  u  the  mean  value  of  u: 
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u(g)  =-^^^/p^  u(x,y)dx. 


Note  that  the  integration  is  over  a  part  of  a  horosphere  centered  at  ~, 
Now  integrate  the  wave  equation  over  F^ : 


^tt  -  y^  "yy  +  ("-2)  y  9y  u  -  (^)^  u=  0 


Introducing  w  =  y^"  '■/i^   ^  and  y  =  e^  as  new  variables  we  obtain 


"tt  -  "ss  =  0- 


We  now  define  the  j*"   component  of  the  translation  representation  as 


^t  -  "s'  i-e- 


nl  n-ls  -X 

R+J{u,u^}  =  c(¥2~^  u^  -  eJg  -2-"^' 


In  part  I  of  [26]  we  show  tliat  R^  is  a  partial  translation 
representation  of  H  for  the  group  Up(t);  in  part  II  we  prove  the 
completeness  of  this  representation.  R,  is  called  the  outgoing 
representation.  One  can  define  quite  analogously  the  incoming 
representation  R._.  The  relation  of  the  two  is  the  scattering  operator 
Sp,  a  notion  introduced  by  Faddeev  and  Pavlov,  [13].  As  pointed  out  in 
Section  A  of  [24],  the  scattering  operator  is  nontrivial  already  for 
the  case  of  F  =  id,  i.e.  for  the  translation  representations  Rg  over 
all  of  H  .   This  Is  in  sh.irp  contrast  to  the  Euclidean  case. 

We  wish  to  emphasize  that  the  translation  representations  are 
constructed  here  in  purely  geometrical  terras,  i.e.  in  terms  of 
Integrals  over  horospheres. 
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3.  The  Zero  Disperslun  Llmtt  for  the  KdV  Equation 


The  equation  in  question  is 


u^  -  6uu^  +  e  2  ^^^^  ^  0  , 


and  the  question  ander  discussion  is  this :  if  the  initial  values   of   u 
are  fixed, 

u(x, 0:  € )  =  u(x)  , 

how  does  the  solution  u(x,t;  e)  behave  as  e  tends  to  0? 

When  we  set  e  =  0  in  the  equation,  we  obtain  the  reduced  equation 

u^  -  6uu^  =  0  . 


This  equation  has  no  solution  for  all  t,  only  in  the  interval 


t^,=  (6  Min  u^(x))  ^      tj  =  (6  Max  u^(x))"^ 
x  x 


here  u(x)  is  the  initial  value  of  u.  It  is  reasonable  to  surmise  that 
for  t  in  (t^,  t^) ,  u(x,t,G)  tends  as  e  ->•  0  to  the  solution  of  the 
reduced  equation.  ^■^at  happens  when  t  lies  outside  this  interval? 
Numerical  experiments  indicate  that  over  some  part  of  the  x-axis, 
u(x,t,e)  is  oscillatory.  As  e  tends  to  0,  the  amplitudes  of  these 
oscillations  remain  finite,  their  wavelengths  are  of  order  e.  Clearly, 
if   we   can  talk  of  a  limiting  behavior  as  e  tends  to  0,  this  limit  can 


20 


exist  only  in  the  weak  sense,  e.g.  in  the  sense  of  distributions. 
This  Indeed  is  the  case;  the  speaker  and  CD.  Levermore  have  shown, 
[22]  ,  that  u(x,t,e)  tends  in  the  sense  of  distribution  to  a  limit  u, 
provided  that  the  initial  value  u(x)  is  nonpositive  and  tends  to  zero 
so  fast  that  /x  u(x)  dx  is  finite.  These  papers  show  not  only  the 
existence  of  a  limit  but  give  a  fairly  explicit  formula  for  the  limit 
u.  For  simplicity  we  take  the  case  when  u(x)  has  a  single  minimum; 
then 


u(x,t)  =  8^^  Q*(x,t)  , 


where 


Q*(x,t)  =  Mn     Q(4';x,t) 
0  <  ij;<4) 

Here  Q  is  a  quadratic  functional  of  u: 


IT  TI 


L  the  linear  integral  operator 


L^(v)  =1/  log  1^1  +(m)  du 


The   functions  admissible  in  the  minimum  problem  are  restricted  to 
lie  between  0  and  (j)  ,  where  (}>  is  defined  in  terms   of   the 
initial  data  as  follows: 


^{r^)=Rej        ^,^%  ..„  0<n 

(lu(y)l-Ti2)l/2 
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Re  /  ^-Jl 

(|u(y)l-u2)l/2 


The  function  a  appearing  in  the  linear  term  depends  linearly  on  x  and 
t,  and  is  defined  as  follows: 

a(n,x,t)  =  nx  -  4n-^t  -  0^(n  ) 

where  0^(n  )  is  a  function  defined  in  terras  of  the  inital  data.  Using 
the  KdV  equation  it  follows  that  also  u^  (x ,  t ;(_)  has  a  limit  u^  -l^^  ^^^ 
distribution  sense,  and  that 


Uj.  =  3ux  . 


Multiplying  the  KdV  equation  by  u  and  rewriting  the  resulting   equation 
as  a  conservation  law  shows  that 


lim  u"  +  -7-  e^  u^  =  u"^ 
e  -f-0      ^ 


exists  in  the  distribution  sense,  and  that 


„2    /   3 
u.  =  4  u^ 


Combining   this   with   the   explicitly   form   the  u  leads  easily  to  the 
formulas 


2      1^       *         '>,  1?* 


The  nvlniimim  problem  dc-fiidr^g  Q*  Is  a  so  cnlU-.d  quadrat.lc 
progranuning  problem;  it  turns  out  that  it  can  be  solved  explicitly.  To 
see  this  it  is  convenient  to  extend  the  functions  ^  admitted  in  the 
nunimum  problem  for  all  real  n  as  odd  functions;  as  a  result  we  may 
replace  the  kernel  of  the  operator  L^  by  log  In-vl-  One  can  show  that 
L  is  negative  definite,  and  that  it  is  related  to  the  Hilbert  transform 
H  as  follows: 


9^  L^H 


We  extend  now  ^i  to  the  upper  half  of  the  complex  n  plane  as  a  harmonic 
function  that  vanishes  at  - ;  ^  can  be  regarded  as  the  real  part  of  an 
analytic  function  of  Hardy  class.  The  variational  condition  for  the 
minimum  problem  can  be  regarded  as  prescribing  the  real  and  imaginary 
parts  of  the  function  on  complementary  subsets  of  the  real  axis;  for 
details  we  refer  the  reader  to  [22].  Suffice  it  here  to  say  that  the 
resulting  formulas  for  u  show 

i)  for  t  in  (t"   t^),  u  is  a  solution  of  the  reduced  equation,  and 

that   in   this   time  interval  the  convergence  of  u  to  G  takes  place  not 

2 
only  in  the  sense  of  distributions  but  for  each  t  in  the  L   sense  m  x. 

ii)  for  t  outside  the  interval  (t",  tj) ,  u  can  be  described  by 
Whitham's  averaged  equations,  or  by  the  more  general  equations  of 
Flaschka,  Forest  and  McLaughlin  based  on  multiphase  averaging. 

iii)  for  t  tending  to  «> ,  u  decays  like  t~  ;  more  precisely 
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a(x,t)  =  -  J_4,d/^)  +  o(t-i) 

Zti    2   t 


for  0  <  x/t  <  Am,  where  m  =  max-  u(x).   Outside  this   range   u(x,t)   is 
0(t-2). 

The  formula  for  u(x,t;£)  is  obtained  from  Gardner,  Greene,  Kruskal 
and  Miura's  solution  of  the  KdV  equation  by  the  scattering  transform. 
VJe  trace  carefully  the  manner  in  which  this  solution  depends  on  e  ,  and 
show  that  as  e  tends  to  zero,  it  has  a  limit  in  the  sense  of 
distributions.  The  nonpositivity  of  the  inital  data  makes  the  GGKM 
solution  of  the  KdV  equation  particularly  simple.  The  more  difficult 
case  when  u  is  nonpositive  has  been  handled  by  Venakides ,  [39]. 
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